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We systematically investigate different versions of variational perturbation theory by forcing not 
only the first or second but also higher derivatives of the approximant with respect to the variational 
parameter to vanish. The choice of the highest derivative version turns out to be the most successful 
one for approximating the ground-state energy of the anharmonic oscillator. It is therefore used to 
determine the critical exponent a of the specific heat in superfluid 4 He in agreement with the value 
" " 1 ", measured in recent space shuttle experiments. 
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(N ■ I. INTRODUCTION 

The perturbative treatment of quantum-statistical or field-theoretical problems renders in general results in the 
form of divergent infinite power series in some coupling constant g. Typically, the coefficients of these series grow 
factorially with alternating signs leading to a zero convergence radius. Various resummation schemes may be applied 
to obtain finite results for all values of the coupling constant g, even in the strong-coupling limit g — > 00 (for an 
overview see Chap. 16 of the book of Kleinert and Schulte-Frohlindc Q and the references therein). Most successful 
I ' is a recent systematic development by Kleinert [|J , extending the variational method of Feynman and Kleinert Q 
which was set up for calculating the effective classical potential in quantum statistics. It has been thoroughly tested 
^\ ■ for the ground-state energy of the anharmonic oscillator and shown to converge exponentially fast and uniform to 
the correct result pj. This was encouraging enough to apply the method also to divergent series which arise from 
renormalizing the </r-theory of critical phenomena fj],||-[| , where the perturbation coefficients are available up to six 
and partly to seven loops in d — 3 [^0|-p^| and up to five loops in d — 4 — e. dimensions |l3| . The method yielded 
finite results with a smooth dependence on the order. Furthermore, the theoretical results are in excellent agreement 
with the only experimental value available so far with an appropriate accuracy, the critical exponent a governing the 
, behaviour of the specific heat near the superfluid phase transition of 4 He @|l5| ■ 

Let us briefly recall the method. Consider some function f(g) which is perturbatively obtained for a small coupling 
£h constant g as the divergent weak-coupling series 
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where the c„ denote the respective expansion coefficients. Kleinert's variational perturbation theory g replaces the 
series ([□]) byQ 



N N-n 



n=o k=a 

By doing so, the newly introduced parameters N , p, q , w are specified as follows. N represents the order of the 
expansion and will be increased step by step as far as the knowledge of the weak-coupling coefficients c n permits. 
The parameters p and q determine the asymptotic behaviour of the function f{g) in the strong-coupling limit g — > 00 
according to 

00 

fig) <? W X><? <*. (1.3) 
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where the bj denote the strong-coupling coefficients. In quantum statistics the parameters p and q are usually known, 
e.g. from dimensional arguments, whereas in statistical field theory they are related to unknown critical exponents, 



1 In contrast to the standard notation |^| the parameter q in Eq. (1.2) has been chosen in such a way that it coincides with 
the critical exponent Q describing the approach to scaling in the field theory of critical phenomena M. 
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so they have to be determined self-consistently from the wea k-coupling coefficients c„ |||[ . Finally, as the variational 
parameter w is introduced such that the approximant ( |l . 2h w ill not depend on it in the limit TV — > oo, it should be 
determined according to the principle of minimal sensitity ]16[. This principle is, however, no clear-cut mathematical 
statement and may therefore be implemented differently, giving rise to different versions of Kleinert's variational 
perturbation theor y. H ere we investigate and compare the versions, which define minimal sensitivity by the vanishing 



of a derivative of (1.2) with respect to w to some order k. Such a version is considered to be good, if its results 
converge well to the true value in the quantum-statistical case, especially when p is known, but q is taken to be 
unknown, since this will anticipate the situation for the field-theoretical application. Another desirable feature can 
be seen in a very smooth dependence on the order TV, since this will greatly enhance the possibility to extrapolate 
field-theoretical results to TV = oo. Finally, the simpler the version's defining prescription, the more easily it may be 
generalized to field-theoretical applications. 

To begin with we focus our attention to the versions of variational perturbation theory in the strong-coupling limit 



oo. Comparing the Eqs. fll.2j) and fll.q) gives, for the leading strong-coupling coefficient bo, the expression 



A' 



to \ N-n J 

where the variational parameter w is still to be optimized. For the ground-state energy of the anharmonic oscillator 
potential 

V(x) = ^x 2 +gx 4 (1.5) 



the leading strong-coupling coefficient has the numerical value 17 



b = 0.667 986 259 155 777 108 270 96... (1.6) 

and dimensional arguments require q = 2/3, p = 1/2. The weak-coupling coefficients c„ are derived in the Appendix 
(see Tab. |v|). 



II. THE HIGHEST-DERIVATIVE VERSION 



Traditionally, the expression (1.4) for bo is made stationary by forcing its first or second derivative with respect to 
the variational parameter w to vanish, depending on whether the order TV is odd or even. Here we investigate whether 
some higher derivatives may be used instead. To this end we consider, for every order TV, all derivatives d k bg /dw k 
with k = 1, . . . , TV + 4 and determine all real po sitiv e zeros w. For each of these zeros, the quantity log \bo — bo\ 
measures the quality of approximating bo in Eq. ( |l.6| ) by b^ in Eq. (1.4) (see Fig. |l|). 

Optimal results are obtained for k — TV. Moreover, the prescription is unique, since for k = TV there exists only 
one real positive zero w of d b$ /dw k . It should be noted that for k = TV — 1 there never exists a real positive zero. 
Another unique and almost optimal prescription would be k — TV — 2 which works well for all TV > 2. The results for 
k = 1 and k = 2 are of comparable quality with the values for k — N. 

Although only one of them has a real positive zero up to TV = 5, the equations for k = 1 and k = 2 tend to have 
more than one solution for larger TV, the one with the largest value for w always producing the best result. But there 
is no simple rule to decide whether k = 1 or k = 2 gives the better value. The initial indication of using k = 1 for 
odd and k = 2 for even orders does not carry through beyond TV = 6. 

In Fig. H we compare the results for k = TV with the optimal branch chosen judiciously from the equations for k = 1 
and k — 2. We also show a quadratic approximation to the k = TV — 2 equation where only the last three terms of 
the condition d 2 bo /dw N ~ 2 = are kept, leaving us with a quadratic equation for the variational parameter w. It 
can be seen that this approximation is of high quality, as was to be expected. 



This may be understood as follows: the coefficients of the w-expansion of 6q in Eq. (1.4) alternate in sign and grow 



faster than n\, which is the same for the weak-coupling coefficients c n from which the former are derived. Therefore 



the series (1.4), or any derivative thereof will preferably become zero, if neighbouring terms in the sum nearly cancel 
each other. The failure of such cancellation among some high order neighbours cannot easily be compensated by 
lower order terms because of the large order behaviour of the coefficients. Hence the vanishing of the derivative for 
k = TV — 2 can well be approximated with little loss of accuracy by choosing w in such a way, that only the three 
highest order terms cancel each other. Thus we conclude that k = N — 2 represents the highest non-trivial derivative 
version for applications in statistical field theory, where we usually have p — 0, so the strong-coupling coefficient b^ 



in Eq. (1.4) is a polynomial in w. It becomes the simpler the higher derivatives are being used. 
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FIG. 1. The approach of bo to its true value &o is measured by the quantity log \ba — b \ which is plotted over k, where k 
appears in the condition d bo /dw k — determining the variational parameter w. Some typical values for the order N have 
been depicted with different symbols T and o used to plot odd and even k, respectively. Optimal results are obtained for k = 1 
or k — 2, but with no simple rule as to which of the two has to be chosen. Taking k = N, however, the same quality is obtained 
with no ambiguity of choice. Also for k — N — 2, the outcome is very reasonable. Note that for some of the lower values of k 
there is more than one solution and that for some k, like e.g. for k = N — 1, there is no solution at all. 
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FIG. 2. The error log 1 60' — bo\ plotted over AT 1 / 3 as obtained by various methods: • for k 
approximation to k = N — 2, and □ for optimal values from k = 1 or k — 2, respectively. 
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FIG. 3. For every N several g-values are obtained by setting the fcth derivative of the strong-coupling limit of 
= gf'(g)/.f(g) with respect to the variational parameter w equal to zero. For each of these self-consistent values of 
q the error log \b$ — bo\ is plotted over N 1/3 . Smaller k give better results, but the highest fc-values up to k = N — 1 are also 
tolerable. 



III. SELF-CONSISTENT DETERMINATION OF THE PARAMETER Q 



In this section we evaluate again the leading coefficient bo in the strong-coupling expansion for the ground-state 
energy of the anharmonic oscillator. But this time the parameter q is not fixed to its proper value q = 2/3. Instead 
we draw this information from the same weak-coupling coefficients c„, whereas p is set to its known value p = 1/2. In 
order to determine q we need an additional equation. Thi s us ually results from a self-consistent reasoning by evaluating 
the logarithmic derivative of the strong-coupling series (F3) of f(g) in the strong-coupling limit g — > oo p,0|: 
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(3.1) 



Defining F( g) :— gf'{g)/f(g), the weak-coupling coefficients 7„ of F(g) are determined from the corresponding c„ of 
f(g) in Eq. (1.1 ). The strong-coupling limit of F(g) can then be constructed with the help of variational perturbation 
theory by using 7„ instead of c n in quite an analogous way as before. Since F(g) tends to the constant value q/2 
in the strong-coupling limit g — > 00, its parameter p has to be zero, whereas q must have the same value as before. 
Thus, from the Eqs. (1.4) and (3.1), we obtain for the strong-coupling limit of F(g) 



n=o \ N -n 




q 

2 



(3.2) 



where we determine the variational parameter w on the left-hand side by demanding that the kth derivative with 
respect to w must vanish, k being chosen appropriately. This set of equations can be rearranged into two polynomials 
of q and w which have to vanish simultaneously. The complete set of complex zeros becomes legion even for moderate 
orders N, their number growing with TV 2 . Constraining ourselves to real and positive solutions for q in the expected 
range 0.5 < q < 0.9, we find, in contrast to the previous case, a very regular behaviour. Exactly one unique solution 
exists for all N up to N = 36 and for all k if N + k is an odd integer, whereas there is no solution within this range if 
N + k is even. The solutions for different k approach the value q = 2/3 closer for smaller k. As k increases, the loss 
of accuracy, however, is tolerable, as can be seen from Fig. [| where again the error log \b$ — b \ is plotted over TV 1 / 3 , 
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the approximation 6q being obtained as before with the highest derivative version, but using the self-consistently 
determined q instead of the exact q — 2/3. Some loss of accuracy is the price to pay, if we want to make one of both 
equations linear by choosing k = N — 1, which is now the highest non-trivial derivative. Notice that the penalty is 
relatively low for small orders N, which in field-theoretical models are the only ones available at present. 



IV. CRITICAL EXPONENT a FOR LIQUID HELIUM 

After having analysed the highest derivative versions of variational perturbation theory by the example of the 
ground-state energy of the anharmonic oscillator, we apply this method to liquid 4 Hc. In particular, we consider 
its superfluid state near the transition point T c , in order to calculate the critical exponent a governing the power 
behaviour \T C — T\~ a of the specific heat. Within the framework of the 4 -theory of critical phenomena, the superfluid 
phase transition of 4 He is described by a complex order parameter field <f) whose bare energy functional is of the 
Ginzburg-Landau type and reads in d = 3 dimensions Q : 

E[<f> B ] = J d 3 x^[dM^} 2 + l<^) + ^9B[^ B (x)] 2] j . (4.1) 

By calculating the Feynman diagrams, one encounters divergencies which are removed by analytic regularization |L8| |. 
A subsequent renormalization leads to renormalized values of mass, coupling constant and field which are related to 
the bare input quantities by the respective renormalization constants Z m , Z g and Z$: 

m\=m 2 Z m Z-\ g B =gZ g Z- 2 , <\> B = ^Z 1 / 2 . (4.2) 

In the literature, one finds expansions for these renormalization constants and for certain logarithmic derivatives, 
the so-called renormalization group functions, up to six and partly up so seven loops |l0|-|l2]]. All these expansions 
can be written in terms of the dimensionless bare coupling constant g B = g B /m. For instance, one obtains, for the 
dimensionless renormalized coupling constant g = g/m, the series 



g(g B ) = J2 K n9 B - (4-3) 

n=0 

The logarithmic derivative of the square mass ratio 

Vm(g B ) = - „ m - B (4.4) 

diogg B 

reads correspondingly 

oo 

Vm(gB) = J2 Xn 9B- (4-5) 

n=0 

Its weak-coupling coefficients K n and A„ have been listed in Tab. [j], respectively. 

In the field-theoretic description of second-order phase transitions, the renormalized square mass m 2 vanishes near 
the transition point T c . Thus the approach to scaling corresponds to the str ong- coupli ng l imit of the dimensionless 
bare coupling constant g B |j]J^-^]. Expecting to obtain finite results for both (4.3) and (4.5), when the dimensionless 



bare coupling constant g B tends to infinity, the strong-coupling expansions for g(g B ) and rj m {g B ) should read 

oo 

g(g B )=g + J2^9 B jn , (4.6) 



3=1 



Vm(gB) ^Vm + ^2ijg B n ■ (4.7) 



5 



TABLE I. Weak-coupling coefficients for gigs) and J7 m (<?s) from 6- and 7-loop perturbation theory ]H)|-[12j (above the 
tripple line) and from extrapolations based on instanton calculations pjlq] (below the tripple line). 



Weak-Coupling Coefficients 




rim = £ A„.gs 


9 = E n n g B 


n 


An 













1 


4 


1 


2 


—46.814 814 814 814 82 


-10 


3 


667.389 693 519 318 6 


120.148 148 148 148 1 


4 


-10 792.618 387 448 09 


— 1 642.256 264 617 284 


5 


191274.332 379 0051 


24 816.045 615 887 86 


6 


—3 644 347.117 315 811 


—407 363.539 559 348 


7 


7 378 080.984 900 002 




7 




7.180 326 000 784143 x E° 


8 


-1.575 312 899 817 985 x E" 


-1.347 981388 551665 x £ 8 


9 


3.529 433 822 947 775 x E w 


2.679 259 494 762 891 x E' J 


10 


-8.269 004 838 427 051 x E 11 


-5.612 792 935 327 37 x E w 


11 


2.020 940 372 700 199 x E VA 


1.234 985 211831956 x E rz 


12 


-5.143 391 961 273 287 x £ 14 


-2.846 297 378 479 211 x E y6 


13 


1.360 628 154 286 311 x E ui 


6.837 009 909 070 767 x E 14 


14 


-3.734 409 500 947 708 x E u 


-1.699 763 366 416 082 x E lb 


15 


1.062 320 774 475 501 x E L " 


4.354 392 675 760 15 x E 1 ' 


16 


—3.132 670 751 194 59 x E 


— 1.151 470 001 1d3 d25 X E 


17 


9.584 565 077 380 565 x E A1 


3.166 449 929 095 059 x E AV 


18 


-3.044 620 939 704 96 x E zi 


-9.135 728 488 698 395 x E' Z1 


19 


1.004119 260 969 699 x E A '° 


2.780 684 275 579 437 x E 1A 


20 


-3.435 379 417 189 77 x E^ 


-8.925 916 735 892 876 x E M 


21 


1.217 685 465 986 233 x E m 


3.006 538 451329 012 x E zh 


22 


-4.465 282 353 357 697 x E™ 


-1.055 6247 984 629 23 x E™ 


23 


1.691830 517 012 908 x E dl 


3.841 496 075 948 586 x E*' 


24 


-6.616 032 677 503 413 x E az 


-1.443 458 944165 208 x E al 



TABLE If. The critical exponent Q as obtained from the strong-coupling limit of the first or the second logarithmic 
derivatives of the weak-coupling series for g(gs) and rj m {gB), respectively. 



SI- Values 


N 






0(?7m, 1) 


n(?? m ,2) 


2 


0.730 495 


0.735 397 


0.715 930 


0.721 303 


3 


0.751 627 


0.751 166 


0.714 270 


0.712 744 


4 


0.757 596 


0.757 762 


0.703 544 


0.700 880 


5 


0.763 865 


0.763 975 


0.705 086 


0.704 576 


6 






0.714 586 


0.714 540 



TABLE III. The critical exponent a for the specific heat of liquid helium, derived from the perturbation expansion for rj, 
using the highest-derivative version of variational perturbation theory. 



a- Values 


N 






2 


0.490 834 


+0.0121 


3 


0.513 786 


-0.0185 


4 


0.522 480 


-0.0304 


5 


0.522 651 


-0.0307 


6 


0.519 592 


-0.0265 



G 



Note that strong-coupling expansions like (4.6) and (4.7) are governed by one and the same critical exponent fl [Ej, 
a common feature of second-order phase transitions. Thus from either of these series fi may be extracted in various 
ways, takin g lo garithmic derivatives with respect to the dimensionless bare coupling constant gs- The strong-coupling 
behaviour (4.6) of the dimensionless renormalized coupling constant .9(35) implies, for instance, 



d log g(g B ) 
hm — = = , 

hm d y<g*) = _ n 



(4.8) 
(4.9) 



with corresponding results for rj rn (g~B)i as can easily been seen from Eq. (4.7). Forcing these relations upon the 
weak-coupling expansions as well, we obtain equations to which the strong-coupling limit of variational perturbation 
theory can be applied. We use its highest derivative version as explained above and find good agreement between the 
fijv-values from each of the four equations (see Tab. |l]). 

In fact they are supposed to approach the same limit as N — > 00. We notice, that the members of each pair obtained 
from either the r\ m - or the g- expansion are almost the same, whereas the different pairs have not yet converged so well 
for N < 6. Even though the data from the g-expansion look more promising because of their smoother behaviour, 
we still use the results fijv from the ^-expansion alone to calculate the critical exponent a, and do not use any 
information from the ^-expansion at this stage. We evaluate the critical exponent r\ m from 



JV 



(-l) iV > X n (-W N 



n=0 




(4.10) 



with the variational parameter wn fixed by 



w N 



X N Nn N 



(4.11) 



Here the Qn are supplied from the fl(r] m , l)-column of Tab. 0. The corresponding results are shown in Tab. [II along 
with the corresponding values for the critical exponent a of the specific heat which follows from 



l-2ry„ 
a = o 

These are values quite close to the experimental data a cxp 
method used. 



(4.12) 



-0.01056 ± 0.00038 |5J, thus giving support for the 
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FIG. 4. The critical exponent £1 derived from the strong-coupling limit of the first and the second logarithmic derivative 
of the weak-coupling expansion of g(gB) (A, V) and of »7m(gs) (■, □), respectively. The solid lines are extrapolations with 
fioo = 0.8 fitted to the g- and ?7 m -points in the reliable interval 3 < N < 12, respectively. 
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V. IMPROVEMENT BY EXTRAPOLATION 



Higher order perturbation coefficients for and r) m (gB) are known approximately (see Tab. |lj). They have 

been derived from instanton calculations, which were fitted to the weak-coupling data Extending the previous 

calculation on this basis up to the order of N = 24, good convergence of all four different sets of the VLm is found. 
We notice, however, a kink at about N = 12 in Fig. [|, which strongly suggests, that the extrapolation is no longer 
reliable beyond this point. For large N the results rijv have the asymptotic form Jl|,|| |j 

M N = ftoo -a exp (-bN 1 -^) , (5.1) 

where the constants a and b are different for the r] m - and the g- expansion. Unfortunately, the fit is extremely 
insensitive to the value of floe', all values in a broad interval around fioc = 0.8 would fit the data up to N — 12 from 
both series with very much the same quality (see Fig. ^) . In Fig. || the resulting a- values are plotted against the 
asymptotic form 

a N = - c exp (-d N 1 ^ 1 ^) , (5.2) 

for ctoo = —0.0106 and (1^ = 0.8. Up to the kink at N = 12, the data seem to approach the experimental number 
excellently. 
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FIG. 5. The critical exponent a of the specific heat of liquid 4 He near the superfluid phase transition plotted against the 
order N. The straight horizontal line represents the experimental value of a exp — —0.0106. The curved line is an extrapolation 
of the results from the range of 3 < N < 12. 

TABLE IV. Comparison of quality for the highest-derivative version and for the traditional method, based on convergence 
of the critical exponent Q obtained in two different ways from two different sets of data. 



an 


data 


method 


2 < N < 6 


7 < N < 12 


13 < N < 24 


9(9 b) 


high 


0.002 159 


0.000 570 


0.017 682 


9 (9 b) 


low 


0.002 245 


0.000 522 


0.000 274 


r?m((/s) 


high 


0.002 749 


0.000 258 


0.016 251 


T]m{gB) 


low 


0.008 214 


0.008121 


0.005 915 


both 


high 


0.046 663 


0.040 836 


0.037 912 


both 


low 


0.046 549 


0.057049 


0.055 477 
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VI. QUALITY CHECK FOR THE HIGHEST DERIVATIVE VERSION 



The four different methods, discussed before to extract the critical exponent £1 from the data, should all give results 
which converge to the same limit as N — > oo. We show now that the highest derivative version reveals a better 
convergence than the traditional method which uses the first or second derivatives instead. The details are reported 

1 from different methods is listed for specific ranges of 



in Tab. IV, where the variance an between the values Sl N 
N: 



\ 



1 u N=N { 



(6.1) 



In the first row of T ab. IV we list the O-values extracted from the weak-coupling expansion ( |4.3| ) of 5(35) according 
to the methods (4.S), (4.9) using the highest derivative version. The same quantity evaluated in the traditional way, 
i.e. using the first derivative for even N and the second derivative for odd N, is shown, for comparison, in the 
second row. It can be seen that both methods are of the same quality in the range of 2 < TV < 6, for which the 
calculation is based on the weak-coupling coefficients from perturbation theory, and also in the range of 7 < N < 12, 
based on extrapolated values. Only for N > 12, where the extrapolation is no more reliable, anyhow, the highest 
derivative version seems to be worse than the one obtained by the traditional method. However, in order to obtain 
these apparently better results using low order derivatives, the rules of the game had to be changed somewhat, taking 
the second derivative throughout for even and odd order N, as soon as N becomes larger than 12. 

The next two rows show the corresponding comparison based on the weak-coupling expansion (|4.5| ) of ?7 m (gs). 
Here it turns out that the highest derivative version gives results of significantly better quality than the traditional 
method, up to N = 12. 



Finally, in the last two rows the devia tion s of the f2-values from the weak-coupling expansions (4.3) of 3(33) and 



(4.5) of rj m (j}B), extracted according to (4.8) are shown in comparison for the highest derivative version (above) and 
the traditional method (below). It can be seen, that both methods are more or less comparable in quality, the highest 
derivative version giving slightly better results. 



VII. APPENDIX 



The perturbative solution of the time-independent Schrdinger equation 

H9 m := (H + gV)V m = E m y m (7.1) 

usually leads to infinite sums for the expansion coefficients of the energy eigenvalues E m (g) and the state vectors ^ m (g)- 
But if the perturbing potential V happens to be band-diagonal in the basis of the eigenvectors of the unperturbed 
Hamiltonian Hq, then all these sums become finite and the eigenvalues E m (g) can be determined recursively to all 
orders. 

Here we assume that the spectrum , eg , e 2 ' , . . . , eg™ 1 ' , ... of the unperturbed Hamiltonian Hq is non-degenerate 
and denote the corresponding eigenvectors by |0), |1), |2), . . . , m), . . . Following the usual procedure, the energy eigen- 
values E m (g) and the state vectors ^ m (g) are expanded in powers of g: 

n,i 

£ m (<?) = E e 5 m V, ( 7 - 3 ) 
3 

where the a n are inserted into the definition, for later computational convenience. They will only show up in the 
coefficients of the state vector fy m , but not in the expansion of the energy eigenvalues E m which we are looking for. 
Without loss of generality the normalisation of the state vectors is chosen such that (^ m \m) = a m to all orders 
so that we have 

7K = ko lW = S k , m ■ (7.4) 



Inserting the ansatz (7.2), (7.3) into the Schrdinger equation (7.1), projecting the result onto the base vector (fc|a& 



and extracting the coefficient of g J , we obtain the relation: 
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with the matrix elements 



n j 



V k , n := (fc|V|n) 



(7.5) 



(7.6) 



For i = this equation is identically satisfied. For i ^ we get the following two relations, one for k — m and the 
other one for k m: 



M _ \ " (m) w 
t i — 2-^i 'm+n,i-l vv m,n 
n 

1 



(m) 
Ik A 



i-1 
3=1 



e j 'k,i—j 



(7.7) 
(7.8) 



In these expressions all sums are finite if the potential V is band-diagonal such that the augmented matrix elements 



^Fm,n • * / m,m+n" 



(7.9) 



are different from zero only for — d < n < d, where d is some finite number. For any quantum number m these relations 



may now be applied recursively with respect to the order, leading to e 1 m ' ) , 7^™ , 7^, ...,£^"',7^ , 

Notice that for given m and r only a finite number of the 7I"!' is non-zero. 

Considering in particular the anharmonic oscillator with the unperturbed Hamiltonian Hq = (p 2 + x 2 )/2 and the 
pert urbing pot ential V = x 4 , we may choose a k = \/2 k jk\ to get the non- vanishing augmented matrix elements from 
(7.6) and (7.G). An easy way to calculate the augmented matrix elements starts from the ones for V — x, where we 
have 



(m) (m) 



. . in turn. 



Wl_ x =k/2, 
wi, = 1 • 



(7.10) 
(7.11) 



From here the augmented matrix elements ■ for any polynomial V = x l can be determined recursively by applying 
the rules of matrix multiplication: 



l-l 
id-i 



(7.12) 



Thus we obtain for V — x 4 : 



W fe % = k(k - \){k - 2)(k - 3)/16 , 
W^ 4 _ a = fc(2fc-l)(fe-l)/4, 

^ 4 = 3(2fc 2 + 2fc + l)/4, 

wf 2 = 2k + 3 , 



W, 



kA 



1 . 



(7.13) 
(7.14) 
(7.15) 
(7.16) 
(7.17) 



Being rational numbers they are suitable for recursive calculations even up to very high orders. In Tab. [v| we have 
listed the first ten coefficients ejj?' for the expansion of the ground-state energy E (g) . 
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TABLE V. The first ten coefficients 4 ' of 
oscillator. 



the weak-coupling expansion for the ground-state energy Eo(g) of the anharmonic 



The Anharmonic Oscillator 



k 


(0) 





1/2 


1 


3/4 


2 


-21/8 


3 


333/16 


4 


-30 885/128 


5 


916 731/256 


6 


-65518401/1024 


7 


2 723 294 673/2 048 


8 


-1030 495 099 053/32 768 


9 


54 626 982 511455/65 536 


10 


-6 417 007 431 590 595/262 144 
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